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Abstract 

We analyze neutrinoless double /3-decay (Oi//3/3-decay) mediated by heavy particles from the 
standpoint of effective field theory. We show how symmetries of the Oz^/3/3-decay quark operators 
arising in a given particle physics model determine the form of the corresponding effective, hadronic 
operators. We classify the latter according to their symmetry transformation properties as well 
as the order at which they appear in a derivative expansion. We apply this framework to several 
particle physics models, including R-parity violating supersymmetry (RPV SUSY) and the left- 
right symmetric model (LRSM) with mixing and a right-handed Major ana neutrino. We show 
that, in general, the pion exchange contributions to Oz//3/3-decay dominate over the short-range 
four-nucleon operators. This confirms previously published RPV SUSY results and allows us to 
derive new constraints on the masses in the LRSM. In particular, we show how a non-zero mixing 
angle £ in the left-right symmetry model produces a new potentially dominant contribution to 
Oi^/3/3-decay that substantially modifies previous limits on the masses of the right-handed neutrino 
and boson stemming from constraints from Oz//3/3-decay and vacuum stability requirements. 
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I. INTRODUCTION 



The study of neutrino less double beta-decay (O^/3/3-decav) is an important topic in par- 
ticle and nuclear physics (for recent reviews, see Refs. (jj, |2, y[ ) . The discovery of neutrino 



oscillations in atmospheric, solar and reactor neutrino experiments proves the existence of a 
non- vanishing neutrino mass ODD . While oscillation experiments provide information on 
mass-squared differences, they cannot by themselves determine the magnitude of the neu- 
trino masses nor determine if neutrinos are Majorana particles. If the neutrino sector of an 
"extended" Standard Model includes massive, Majorana neutrinos, then O^/3/3-decay pro- 
vides direct information on the Majorana masses. Indeed, since Majorana neutrinos violate 
lepton number (L) , Feynman graphs such as the one depicted in Fig. ^ are non- vanishing. 
In particular, if the e, Li, t neutrinos have non-vanishing Majorana masses, an analysis of 
OvfiP coupled with data from neutrino oscillations provides limits on the absolute value of 
these light neutrino masses 

Neutrinoless /3/3-decay can also be a probe for heavy mass scales. For example, in the 
left-right symmetric model [3,0,0], a heavy right-handed neutrino also contributes to the 
process; it can even be dominant depending on the values of the elements of the mixing 
matrix. Thus, 0^/3/3 can be a tool for the exploration of energy scales beyond the electroweak 
symmetry breaking scale. Alternatively, the L-violating interactions responsible for Qv/3/3- 
decay may not involve Majorana neutrinos directly. For example, semileptonic, R parity- 
violating (RPV) supersymmetric (SUSY) interactions, involving exchange of charged-lepton 
superpartners (an example of which is given in Fig. rather than Majorana neutrinos, can 



FIG. 1: a) Qvj3(3 through the exchange of a Majorana neutrino, b) 0vf3f3 through the exchange of two 
selectrons and a neutralino in RPV SUSY. 
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give rise to Oz//3/3-decay |10L IllL Il2j |. Here again 0z//3/3-decay provides a probe of the heavy 



SUSY mass scale and imposes constraints on RPV SUSY parameters 



3- 



Furthermore, 



these alternative scenarios for O^/3/3-decay are relevant for the study of Majorana neutrinos 



since any 0z//3/?-decay mechanism will generate Majorana masses for the neutrinos 

The left-right symmetric model and RPV SUSY are but two of a number of models that 
involve a heavy mass scale App that characterizes the heavy, L-violating physics. Although 
the effects of these mechanisms will typically be suppressed by some inverse power of A^, 
Oz//3/3-decay mediated by light neutrinos can also be suppressed since the amplitude is pro- 
portional to the neutrino effective mass. Thus, it is important to analyze systematically the 
potentially comparable contributions stemming from L-violating mechanisms mediated by 
heavy particles. Since App is far heavier than any hadronic scale that would enter the prob- 
lem, there exists a clear separation of scales in this case. For the analysis of such situations, 
effective field theory (EFT) is the tool of choice. 

In what follows, we systematically organize the Oz//3/3-decay problem using EFT, fo- 
cusingon L-violation mediated by heavy physics (for other efforts along these lines, see 
Refs. |15I.I1(tLI17I|). Since the particle physics dynamics of this heavy physics occur primarily 
at short- distance, one may "integrate out" the heavy degrees of freedom, leaving an effective 
theory of quarks and leptons; these quark-lepton operators in turn generate hadron-lepton 
operators that have the same transformation properties under various symmetries. In this 
work, only the lightest quarks are considered, with the relevant symmetries being parity 
and strong SU(2)lxSU(2)r [chiral SU(2)]. The effective hadron-lepton Lagrangian for this 
theory, C^pp, contains an infinite tower of non-renormalizable operators, which may be 
systematically classified in powers of p/A H , p/ App and An/ App. Here, p denotes any small 
quantity, such as m n or the energy of the dilepton pair and Ah ~ 1 GeV is a hadronic mass 
scale. While the coefficients of the effective operators in are unknown 1 , the symme- 

try properties of the underlying short- distance physics may require that certain operator 
coefficients vanish. 

These symmetry properties can have significant consequences for the size of O^/3/3-decay 
nuclear matrix elements and, thus, for the short- distance mass scale deduced from exper- 



1 The computation of these coefficients from the underlying quark-lepton interaction introduce some degree 
of uncertainty - a problem we will not address in this work. 
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imental limits. Specifically, the hadronic vertices appearing in C^pp will be of the type 
NNNNee, NNiree and iriree, etc. They stem from quark-lepton operators having different 
transformation properties under parity and chiral SU(2); as such, they will contribute to 
different orders in the p/Ah expansion. 

Traditionally, the short-range NNNNee contribution to O^/3/3-decay has been analyzed 
using a form-factor approach |l8| where the finite size of the nucleon is taken into account 
with the use of a dipole form-factor. The form-factor overcomes the short-range repulsive 
core in NN interactions that would otherwise prevent the nucleons from ever getting close 
enough to exchange the heavy particles that mediate Oi//3/3-decay. The disadvantage of a 
form-factor model is that the error introduced by the modeling cannot be estimated system- 
atically in contrast to the EFT approach. A discussion of the NNNNee vertex within the 
framework of EFT will appear later in this paper. 

In contrast to the short range contribution to O^/5/3-decay, the long range contributions 
involve the exchange of pions pj| through the NNiree and irnee vertices. Although these 
long range contributions have been analyzed in the form-factor approach j^j], they are 
more systematically analyzed within the context of EFT because of the separation of scales: 
< Ah <C App. As noted in ref. for example, the matrix elements associated with the 
long range pionic effects allowed under RPV SUSY scenarios can be dominant. However, we 
show that the dominance of pion exchange in Ozz/3/3-decay mediated by heavy physics is a 
more general result not limited to RPV SUSY. These pionic effects can be considerably larger 
than those obtained using the conventional form factor model for the short-range NNNNee 
process. For these reasons, the analysis of the long range contributions to Oz//3/3-decay in 
EFT will be the main focus of this paper. 

The various types of L- violating operators that contribute to the long range contributions 
of O^/3/3-decay appear at different orders in the p/Ah expansion with p ~ m n , and the order 
at which they appear depends on their symmetry properties. It is therefore important to 
delineate clearly the symmetry properties of Cppp for various types of L- violating operators 
and use these symmetries to relate the hadron-lepton operators to the underlying quark- 
lepton operators. Carrying out this classification constitutes the first component of this 
study. In doing so, we also comment on the standard approach to deriving Oz//5/3-decay 
nuclear operators and correct some errors appearing in the literature. 

The second step in our treatment involves deriving 0v(3(3-dec&y nuclear operators from 
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^EFF expressing the rate in terms of corresponding nuclear matrix elements. For any 
/9/9-decay mode to occur, the final nucleus must be more bound than any other prospective 
single /3-decay daughter nucleus. Such /3-forbidden but /3/3-allowed nuclei only occur for 
sufficiently heavy nuclei. Thus, the extraction of the short-distance physics that gives rise 
to Ou (3(3- decay (at present, only upper limits on the decay rates exist) depends on a proper 
treatment of the many-body nuclear physics. Having in hand the appropriate set of nuclear 
operators (for a given L- violation scenario), one could in principle compute the relevant 
nuclear matrix elements. Unfortunately, it is not yet possible to do so in a manner fully 
consistent with EFT. This problem has been studied extensively in the case of the NN 
and three-nucleon systems, where the state-of-the art involves use of chiral symmetry to 
organize (and renormalize) the relevant nuclear operators 
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Out of necessity, 



we follow the same philosophy here. Nonetheless, the organization of various 0z//3/3-decay 
operators based on symmetry considerations and EFT power counting should represent an 
improvement over present treatments of the nuclear problem. 

As a final step, we relate the various nuclear operators obtained from t° different 

particle physics models for L-violation. Doing so allows us to determine which nuclear 
mechanisms dominate the rate for a given particle physics model. For example, in both the 
RPV SUSY and the left-right symmetric model with mixing of the gauge bosons, the rnree 
contribution to the Oi^/9/3-decay amplitude is significantly larger than that of the short range 
NNNNee contribution. In contrast, for left-right symmetric models with no mixing, these 
contributions are of a similar magnitude. We also show how this large iriree contribution 
to Ov (3 /9-decay substantially affects the relationship between the masses of the right-handed 
neutrino and gauge boson including a new correlation between the minimum mass of the 
right-handed neutrino and the Wl — Wr mixing angle. In short, the sensitivity of the 
0i//3/?-decay searches is strongly affected by the symmetry transformation properties of the 
operators contained in a given particle physics model. 

The remainder of our paper is organized as follows. In Section[nJ we classify the operators 
in C^(pp according to their symmetry properties and p/A counting and we tabulate the 
various quark-lepton operators according to the hadron lepton operators they can generate. 
In Section IIIII we use the leading operators to derive non-relativistic nuclear operators and 
compare their structure with those appearing in conventional treatments. In section IIVI 
we work out the particle physics implications under various scenarios, namely RPV SUSY 
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and the left-right symmetric model and compare them to each other. We summarize our 
conclusions in Section fVl 



II. EFFECTIVE Oz//?/?-DECAY OPERATORS 



The classification of the operators in £eff renes on two elements: 

1. The use of symmetry to relate effective lepton-hadron GV/3/?-decay operators to those 
involving quarks and leptons. The relevant symmetries are parity and chiral SU(2). 
Indeed, because the lepton-hadron effective operators are generated from the quark- 
lepton operators through strong interactions, they should retain the same parity and 
chiral structure. 

2. The organization of these effective lepton-hadron operators in an expansion in powers 
of a small momentum p. 

To organize the non-standard model (NSM) operators in powers of p, consider first the 
long range 7r-exchange contributions to GV/3/3-decay of Figs. |2^,b, and c. The fact that pions 



are Goldstone bosons allows us to use chiral perturbation theory |26|,|27[ to classify the NSM 
hadronic operators in terms of a p/An expansion, with Ah = ~ 1 GeV and p ~ m n 
where /„- ~ 92.4 MeV is the pion decay constant. The leading order (LO) quark operators 
should therefore induce effective hadronic operators that do not involve derivatives of the 
pion fields or pion mass insertions 2 , the next-to- leading order (NLO) operators would involve 
a single derivative of the pion field, the next-to-next-to-leading order (NNLO) would involve 



FIG. 2: Diagrams that contribute to Qvj3j3 at tree level. The exchange diagrams are not included. 
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2 At tree level, the pion mass insertions always have the form m 2 . and therefore do not contribute at LO or 
NLO. 
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two derivatives or pion mass insertions and so-on. This approach to Oz//3/3-decay is similar to 
the application of effective field theory to purely hadronic AS* = parity-violating operators 
that was done in and the same notation will be used. 

The power counting for the long-range Oz//3/3-decay operators will involve the chiral order 
of the standard model (SM) operators as well as the chiral order of the NSM operators. For 
the SM operators, these counting rules are as follows: 

• a pion propagator is 0(l/p 2 ) while 

• each derivative of the pion field and the LO strong ttNN vertex are 0(p). 

As for the short range operators (Fig. |21i, the hadronic part is constructed from a 4- 
nucleon vertex. This vertex can also be expanded in powers of the nucleon's 3-momentum. 
However, the chiral counting suggests that the leading O(p ) four-nucleon vertex is already 
strongly suppressed relative to the long range 0z//3/?-decay operators such that the 4-nucleon 
vertex can be neglected to lowest order. Indeed, with these rules, the chiral counting of the 
Oi^/3/9-decay operators of Figs. H^-d are 

Fig. Hi ~ K wn p~ 2 , Fig.Eb,c~iW^\ Fig.Hl- K NNNN p°, (1) 

where the K{ denote the order of the NSM hadronic vertices. In general, the LO vertex 
in each diagram is O(p ), though in certain cases symmetry considerations require that 
the leading order vertex vanish (see below). Thus, the long range O^/3/3-decay operators of 
Figs. Ek, and |2)d,c are enhanced by 1/p 2 and 1/p, respectively, relative to the short-range 
operator of Fig. Eli. In what follows, we will consider contributions generated by all of 
the diagrams in Fig. |21 Since the LO contribution from Fig. |21i is O(p ), we must include 
contributions from Fig. |2^-c through this order as well. Consequently, we consider all terms 
in K nn and Knntt to 0{p 2 ) and 0(p), respectively. 

A. Quark-Lepton Lagrangian 

In order to construct the hadron-lepton operators, we begin by writing down the quark- 
lepton Lagrangian for O^/3/3-decay. This is done by considering all the non- vanishing, in- 
equivalent, lowest-dimension quark-lepton operators that are Lorentz-invariant and change 
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lepton number by two units, 



2 



A/3/3 

+ (o 6 0+ + + + o 7 2 + + + 08 O+ + + o 9 3 + + + + O i O 3 + - + ) e 7 V 

+ (onCi +,/i + o 12 Ot^ + o 13 5 + + +l/i + oi4C 5 + _ + ' M ) e 7/ , 7 V + h.c.} , (2) 



where 



Cft = (&TV<fc)(&A,*«R), (3) 

Cft = (gRr a gL)(gRrV)±(gLr a gR)(g L rV), (4) 

(9 3 a i = (qLT a rqL)(qLr b ^q L )±(q R T a rqR)(qRT%qR), (5) 

4 a ± M = (9ltV<7l =F ? R tV?r)(?ltV - g R r 6 g L ), (6) 

Off = (g L rV?L ± <?rtV<?r)(<?lAr + g R r fe g L ). (7) 

The 9l,r = (w, gOl,r are left-handed and right-handed isospinors and the r's are Pauli matri- 
ces in isospace. When a = b, the operators with subscript +(-) are even (odd) eigenstates of 
parity as can be verified by noting that the parity operator simply interchanges left-handed 
spinors with right-handed spinors. This list of nine operators was arrived at by inspec- 
tion 3 . Other operators that could have been written down are either equivalent to those 
in Eqs. Q to (JJJ) or vanish as shown in appendix [X] in particular, all operators propor- 
tional to ea^e , e , y 5 a fiu e c and e^e c vanish since these leptonic currents are identically zero 
as can be verified with the use of Fierz transformations. Some of these vanishing leptonic 

n 

currents were erroneously taken as non-zero in Ref. [17J|. Similarly, a quark operator, like 
qa ,lu T ± qqa f j il/ T ± q, can be re-expressed in terms of 0^± by applying a Fierz transformation 
despite the color indices since the hadronic matrix elements of four-quark operators only 
select their color singlet part 4 . 

Recalling that fermion fields have mass dimension 3/2, note that the operators appearing 
in Cfogp have mass dimension nine. Therefore, the overall coefficients have dimensions 
[Mass] -5 . In Eq. (J2J), this scale factor is expressed as Gp/App where A/3/3 remains to be 



3 In writing down the Eqs. 1- !I7I) . we suppressed the color indices since EFT only relates color-singlet quark 
operators to hadronic operators. 

4 The projection onto color singlet states introduces a new factor that can ultimately be absorbed in the 

OiS. 
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determined. Derivative quark operators are suppressed by extra powers of A^g and need not 
be considered further. 

The operators in C q Qu pp can be generated by various particle physics models, but not all 
of them are necessarily generated in a single model. For example, the left-right symmetric 
model always involves the product of left-handed and/or right-handed currents, while only 
Of+ and 0%± are of that form. Thus, 0%±, O ^ and Q§± cannot appear in the left-right 
symmetric model. Another example is a minimal extension of the standard model with only 
left-handed currents and Majorana neutrinos; in this scenario, only Of ± could appear. On 
the other hand, these operators all appear in RPV SUSY. This observation will allow a 
classification of these particle physics models later in this paper. 

Since O^/3/3-decay always requires a — b — ±, the O's have definite transformation 
properties . Using the quark field transformation properties under chiral SU(2), 

under SU(2) L xSU(2) R : q L - Lq L , q R - Rq R , (8) 

where the L and R transformation matrices have the form oxp{Pl^r9l,r} and 

e L ,R = \r-e L , R , P L , R = l -{lTl 5 ), (9) 

we derive the transformation properties of the 0^9^ under chiral SU(2), 



/nab 
^1+ " 


+ (q L L^T a rLq L )(q R RW\Rq R ), 


(10) 


/r\ab 
U 2± ~ 


+ (q R R^r a Lq L )(q R R^r b Lq L ) ± (q L L^r a Rq R )(q L L^r b Rq R ), 


(11) 


/r\ab 
w 3± 


+ (q L L^T a YLq L )(q L L^T b ^Lq L ) ± (q R R) T a ^^ Rq R ) (q R R^T b '-f f j i Rq R ) , 


(12) 


/nab,fi 
U 4± ~ 


■* (q L LW a ^Lq L T q R R^T a ^Rq R ){q h L\ h Rq R - q R RW b Lq L ), 


(13) 


/r\ab,fi 
U 5± 


+ (q L L^ a YLq L ± q R R^ r a, -f^ Rq R ) [q\ J L^T b Rq R + q R Rh b Lq L ). 


(14) 



We observe that O^. belongs to the (3^,3^) representation of SU(2) L xSU(2) R (from here 
on, the subscripts L, R are dropped) in the sense that the first superscript a transforms like 
a triplet under SU(2) L while the second superscript b transforms like a triplet under SU(2) R . 
Note that only Of+ belongs to a representation of chiral SU(2). The other s are 

superpositions of operators that have different transformation properties under chiral SU(2). 
This is not surprising since the generators of chiral SU(2) do not commute with the parity 
operator as they involve 7 5 . For instance, (?LT ± 7 /i Q'L)(?L' r± 7^Q , L) changes isospin by two 
units and is a singlet under SU(2) R such that it belongs to (5,1) while (qRT ±r y fJ 'q R )(q R T ± 'j^q R ) 
belongs to (1,5). Hence, Of± belongs to (5, 1) © (1,5). 
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B. Hadron-Lepton Lagrangian 



Let us now turn to the derivation of the iriree vertex from the quark operators. This will 
be followed by a similar analysis for the NNnee and NNNNee vertices. 



1. irnee vertex. 

To derive the hadronic vertex, first consider parity. The product of two pion fields being 
even under parity, only positive parity operators can contribute. Secondly, note that O^' 11 
and 054^ must give rise to an operator of the form 

7r+^7r + e7 M 7 5 e c + h.c. (15) 

A partial integration shows that this operator is suppressed by one power of the electron 
mass, and is therefore negligible. 

Thus, the only terms in C^ v ^ that contribute are: 

P- {Of+e( 0l + o 67 5 )e c + 0£+e(o 2 + o r7 5 )e c + 0^e(o, + o 9l 5 )e c + h.c.} . (16) 

The hadronic operators that stem from these quark operators must have the same trans- 
formation properties and can be written down by introducing the following fields 

x* = = ev a e, x a = (n) 



£ = exp^yr/X) = exp 



i ( _L _L . _ _ 1 



3^0 



r 7i + r Ti H —t 71 



V2f« V V2 



(18) 



vr ± = ^(tt 1 t «vr 2 ), N : Nucleon field. (19) 
v2 

The transformation property of the above fields under parity are 

7i — > — 7r, £^£ f , X£^X£, X a ^X ta , N->-f°N, (20) 

while under SU(2) L xSU(2) R they transform as 

£ - L£rf = U£tf (21) 

X a -> U£ R ] r a L £U ] (22) 

XI -> Utf L\ a L £U ] (23) 

X£ -> tff i?V a i? ^ (24) 

N ^ UN . (25) 



The transformation matrix U only depends on the r's and the pion field. 

At LO (no derivatives), the two-pion operator stemming from the 1+ operator is 

Oft - tr[$g=] = tr[X±Xi + x£x£] = Att V + ■■■, (26) 

J TT 

while the one generated by Of + is 

Of + tr[$g=] = tr[X ± X ± + Xt±Xt±] = -- I^tt* + ■■■ . (27) 

J 7T 

Here, $ x 2± are defined 

*g = *L *R ± *R*L , 

= X ± X ± ± X ]± X^ , (28) 

and the ± subscript refers to the transformation properties of the 's under parity. 

Note that when the traces of <f>f+ and $3+ are expanded up to two powers of the pion 
field, they are physically indistinguishable since the relative minus sign can be absorbed in 
a operator coefficient referred to as a low energy constant (LEC) . 

Now consider the case of the two-pion operator generated by Of+; to LO the hadronic 
operator should be: 

tr [X+X+ + X+X+] = . (29) 

Thus, there exists no (5, 1) © (1, 5) hadronic operator with no derivatives. 
The LO Lagrangian for the iriree vertex is therefore 

q*ee = 3_ | tr [$++]e( a + 6 7 5 ) e c + tr[$r/+ Je^a + b-f 5 )e 

A/3/3 

+tr[$ 2 \ + ]e(a / + 6 / 7 5 )e c + tr[$ 2 --]e c (a / + 6 / 7 5 )e} , (30) 

where a, b, a', b' are LEC's. Note that although there are nominally four LEC's, once the 
traces of the $f's are expanded, there are in practice only two: a — a' and b — b' . 

In contrast to the Oj's, the a,b,a',b' are dimensionful. It is useful to express them in 
terms of dimensionless parameters (denoted in this work by Greek letters) with the aid of 
a scaling rule. In a scaling rule, the hadronic operators are divided by the relevant scales 
such that their coefficients are dimensionless and of a "natural" size. We follow the naive 
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dimensional analysis (NDA) scaling rules given in Ref. 29| and modified here to account for 
the lepton bilinears 5 : 

(SH0(0''(fH 

Justification for this scaling rule is given in Appendix [Bj Note that the scaling factor 
(tt/ ' f^) m is already properly accounted for in the definition of £ and need not be applied 
again in Eq. (JHOj) after expanding the $'s to two pions. For the non-derivative imee vertex, 
we have (k, /, m) = (0, 0, 2) and 

™ = GfAjLg l^-^-g^ + ^ )(f + 7r+7r+e c (A - (3 2 ^)e) . (32) 

Consider now the higher order contributions to the mree vertex. As discussed below 
Eq. (j!5|) . there is no NLO contribution. Hence, Cj^ ee = 0. 

At NNLO, not only do Of+ and Of+ generate two-derivative hadronic operator, but so 
does Of± 

of} -> l -tx [v»x±v,xt + v»x£v„x£] , (33) 

where the chiral covariant derivative is given by 

V^d.-iV^ V M = if(^ t + ^e)- (34) 

The operator "Z^At^r has the same transformation properties under chiral SU(2) as Al,r. 

The only other contribution stems from quark mass insertions that always generate 
squared pion mass insertions. Writing the NNLO contributions directly in terms of pion 
fields, we obtain 

= ^ {d^-d^-e(f3 3 + &7V + mln-ir-e({3 5 + /? 67 5 )e c + h.c.} . (35) 

A/3/3 

Note that the (3*,$ terms constitute corrections to — > /?i,2 + "^105,6 that can be ignored 
in particle physics models where the LO operators contribute since /3i 2 must be measured 6 . 



5 We neglect electromagnetic effects. 

6 As discussed in Ref. |30j , EFT relates the two-derivative mree operator to the 27-plet K ^ 2n decays 
indicating the possible existence of an extra suppression factor beyond that deduced from power counting. 
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2. NNnee vertex 



We analyze the NNiree vertex of Figs. Eh> and Et using similar logic as in the foregoing 
discussion. The LO Lorentz-scalar NNn operator is Nt ± tt t N which is odd under parity. 
Therefore, Of+, and Of+ cannot contribute since they are parity even. As for Of_, 

notice that as in the inree case, the LO contribution (X^X^ — X^X^) vanishes. 

The operator Nr^ir^N can only be induced by Of^. The result is 

C>±± -> N®f±N. (36) 

It is straightforward to verify that X<J? 2 _ X transforms precisely like O^t under 
SU(2) L xSU(2) R . 

In addition, O^'^ and Of± ,fl also generate LO contributions to the NNn operator, 

0^,0^ -> NYl 5 $i-N, 

0^,0^ -> N-f^f-N, (37) 

where 

$f_ ± = (X L ± + X±)(X ± -X ± t), (38) 

as can be checked explicitly by considering the transformation properties under chiral SU(2) 
and parity. The NNnee LO Lagrangian can now be written down, 

c NN,ee = G? [NQ++ N g( c + + ^-f{f X + M^^^V + h.C. } 

A/3/3 

= {Xr + 7r^Xe(Ci + C 2 7 5 )e c + XytCs + C47 5 )r + ^~A^ 7 5 e c + h.c.} ,(39) 

where the Q are dimensionless LEC's introduced using Eq. (|3*T|) with (k, /, m) = (1, 0, 1) and 
where we have expanded the $'s to one pion. 

At NLO, Of±, 0^ and Og= contribute to the XXvr operator, 

Of+ -> X 7 5 $±±X, (40) 

- X^^X, (41) 

±± _> X {y [X±(-zP M X±) - X±(zP M X±)] } X, (42) 

03 d ? ^ # {7 M 7 5 [X±(-zP M X±) - X±(z2^Xr)] } N (43) 
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The first thing to note is that a term like N^irN is subleading because in the non-relativistic 
reduction, the 7 5 couples small and large components of the nucleon spinors. Secondly, we 
observe that Eqs. ()4(J|). (fTTj) and (|4*5)l are physically indistinguishable on shell when expanded 
to one pion and to the order we are considering, as seen from the equations of motion. 
Thirdly, Eq. (|42j) is negligible even at NLO because the equations of motion can be used 
to show that N $ttN is proportional to the electron momentum. Therefore, Of^ does not 
contribute to the NNiree vertex. 

Other contributions to 0(p) include terms normally neglected at LO in the non-relativistic 
reduction of Eq. (jSHJ), namely the terms proportional to £3 and £4 with /j, — 1,2,3 and /1 = 
respectively, where LO and NLO components of the nucleon spinors are coupled. These 
are the only contributions to the NNiree vertex since the m 2 insertions are of 0(p 2 ) and 
excluded as discussed below Eq. (JTJ). Hence, the only new contributions to 0{p) is, 

A/3/3 

where the scaling rule in Eq. (|31J1 was used with (k,l,m) = (1, 1, 1). C^Nnee is subleading 
because the 7 5 couples the large and small components of the nucleon spinors and the result 
is proportional to p/M where M is the nucleon mass and p is the magnitude of the nucleon 
three-momentum. 

3. NNNNee vertex 

To identify the quark operators that generate the O^/3/3-decay four- nucleon operators, we 
insert the hadronic fields X^, X ±:k , X^ ±± in all possible ways into NT N NT' N and use 
their transformation properties under chiral SU(2) to relate them to the Of± . The four- 
nucleon operators are then obtained by expanding these hadronic fields to LO and ignoring 
all contributions from pion loops. Thus, it is not necessary to insert these hadronic fields 
in all possible ways; we only need to show that a particular quark operator can generate a 
particular nucleon operator with the same transformation properties under parity and chiral 
SU(2). 

For example, the LO operator (Nt^N) 2 can be generated by Of+. The latter transforms 
the same way under parity and chiral SU(2) as the hadronic operator 

(NX±N)(NXiN) . (45) 
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TABLE I: Cross-reference table between nucleon and quark operators. The X indicates that the 
quark operator cannot generate the corresponding nucleon operator while the \J indicates that it 



can. 



NNNN ops. 


u i+ 




ot± 


^3+ 




n ±±,/J. 










v 7 


V 


X 


V 


X 


X 


X 


X 


X 


Jl 2+ 


v 7 


V 


X 


V 


X 


X 


X 


X 


X 


Jl 3+ 


V 


V 


X 


V 


X 


X 


X 


X 


X 




X 


X 


X 


X 


X 


V 


X 


V 


X 




X 


X 


X 


X 


X 


X 


V 


X 


V 



At zero pion order, the and both become r 1 * 1 , so that the operator in Eq. (|45jl 
just becomes (Nr^N) 2 . In a similar fashion, it can be easily shown that the following five 
operators 



<tt±± = (N^N) 2 , W±± = (iVr ± 7 /1 iV)(iVr ± 7 /1 iV), ^ = ( Nt^j" N) ( N^j 5 ^ N) 

l 4+ 



fj\±± * = ( A r r ± 7 M A^) (Nt^N) , W^'" = {Nt^^N^N^N), (46) 



exhaust the list of possible LO four-nucleon operators 7 that can be generated by the checked 
Of^ ], s in Table □ 

The LO four-nucleon Lagrangian is therefore given by 



r NNNNee 
^0 



<_Tp 



{ {ii^it + 6*2+ + 6^ 3 + +) ee c + {^nft + + 6*3+ ) n 5 e c 

+ (£7*4+^ + ZsKt-'* 1 ) e 7 5 7 M e c + h.c. } , (47) 



where the £j's are dimensionless. 

In concluding this section, we discuss a few issues that will require future work. The 
first involves the application of EFT to heavy nuclei. As pointed out earlier, no fully consis- 
tent treatment for such situations has yet been developed. In principle, one could imagine 
following a program similar in spirit to the EFT treatment of few-body systems. In that 
case, there has been recent progress in developing a consistent power counting for EFT with 



explicit pions 



251 ] . The approach involves including the LO 7r-exchange contribution to 



7 Since N*y 5 N and (Nj 5 'y ll N)(Nj l j j N) are proportional to p/M, they are sub-leading in the non-rclativistic 
limit. 
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the NN potential, expanding it about the chiral limit (m^ — > 0), and obtaining two-body 
wavefunctions by solving the Schrodinger equation with the chirally-expanded potential. To 
be consistent, operators would also be expanded to the same chiral order as the potential 
and matrix elements computed using the corresponding wavefunctions. This approach ap- 
pears to reproduce the consistent momentum power counting obtained with perturbative 
pions in the 1 S'o channel and the convergence obtained with non-perturbative pions in the 
3 S\- 3 Di channel. In going to more complex nuclei, one might explore a marriage of the chiral 
expansion with traditional many-body techniques (e.g., shell model or RPA), in which case 
one would require a corresponding chiral counting of nuclear operators. In organizing the 
0u j3 /3-decay hadronic operators according to both the derivative and chiral expansion, we 
have taken one step in this direction. For the moment, however, we will have to content 
ourselves with using these operators along with wavefunctions obtained from traditional 
many-body techniques. 

A second issue is the presence of higher partial waves in the two-body transition matrix 
elements appearing in 0z//3/?-decay. A fully consistent treatment would, therefore, require 
that one include the corresponding higher-order operators - a task that is clearly impractical 
at present. Fortunately, in our case, there is reason to believe our qualitative conclusions 
about the dominance of long-range, pion-exchange operators are fairly insensitive to this 
issue. For the cases where the LO iriree are not forbidden by the symmetries of the quark- 
lepton operators, the LO 7r-exchange operator arising from Fig. |2K will always give the LO 
contribution to the transition matrix element, regardless of the partial wave decomposition 
of the two-nucleon initial and final states. In general, then, we expect that matrix elements 
of these operators should always be enhanced relative to those involving the four-nucleon 
contact operators or 7r-exchange operators obtained with higher-order pionic vertices. In- 
deed, some evidence to this effect is given by the computation of Ref . j^] , where the relative 
importance of the LO 7r-exchange operators and short-range operators were compared for 
RPV SUSY 8 . 

Finally, when NNLO and NLO interactions are included at tree level, loop graphs must 
also be included to be consistent with the power counting (examples of which are given 
in Fig. EJ). These loop graphs are handled according to the chiral perturbation theory 

8 However, in that work, the traditional, form factor approach was used to compute short-range effects. 



15 



prescription by which the divergences renormalize the LEC's that multiply the and two- 
derivative iriree vertex of Eq. (|35|). In this context, loop graphs that renormalize the NNiree 
vertex are N 3 LO and can be ignored. Indeed, this can be demonstrated using power counting 



where each loop involves a factor of p 4 while nucleon propagators count as p~ l \3uL 1371 138|. 

When loops are included, new lepton- violating tree level vertices can contribute inside the 
loop graphs, such as the iriririree vertex of Fig. Eh- Other new vertices that could potentially 
contribute at the one loop level are NN iriree and iririree vertices. In short, large number of 
Feynman diagrams may need to be calculated at NNLO. We defer a discussion of such loop 
contributions to a subsequent study. 



To summarize the conclusions of the analysis, Table |H] lists the quark-lepton operators 
that contribute to the various hadron operators at LO. One important result indicated in the 
table is the fact that if the short- distance physics responsible for Oz//3/3-decay belongs to a 
representation of SU(2)lxSU(2)#, only operators that belong to the (3,3) and (5, 1) © (1, 5) 
can generate Oz//3/3-decay and therefore, only Of+ and Of+ can contribute. For example, 
the left-right symmetric model with mixing between left- and right-handed gauge bosons 
induces operators belonging to the (3,3) as well as the (5, 1) © (1, 5). From Table ITU the LO 
O^/3/3-decay operator that contributes in this case is generated by Fig. [2^ and is 0(p~ 2 ). 

Alternatively, consider a short- distance model involving products of two left-handed cur- 
rents or two right-handed currents only. Such a situation arises, for instance, in the left-right 
symmetric model when the Wl and Wr bosons do not mix. For this scenario, only Of+ 
contributes, and there are no LO contributions to the iriree and NNiree vertices. The first 

FIG. 3: a) Example of a graph that renormalizes the LEC's that multiplies the rrv\ and two-derivative 
iriree vertex, h) Example of a new vertex (iriririree) that contributes to §v(3(3 at NNLO. 
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TABLE II: Leading order O^/3/3-decay hadronic-lepton operators generated by the various quark- 
lepton operators. 



0^,3,3-decay ops. 


oft 


±± 




of? 






±±,„ 


±±,„ 


±±,„ 


rnree LO 


V 


V 


X 


X 


X 


X 


X 


X 


X 


7T7ree NNLO 




V 


X 




X 


X 


X 


X 


X 


NNiree LO 


X 


X 


V 


X 


X 


V 


V 


V 


V 


NNiree NLO 


X 


V 


X 




X 


V 


V 


V 


V 


NNNNee LO 


V 


V 


X 


V 


X 


V 


V 


V 


V 



non-zero contributions to the hadronic part of these vertices are given by Eqs. (J3*3*j) and 
()43|) as well as contributions that include m 2 insertions. The resulting contribution to the 
amplitude is O{p ). In this case, both the long- and short-range nuclear operators occur at 
the same order. 



III. NUCLEAR OPERATORS TO LO AND NLO 

In the calculation of the Oz^/3/3-decay amplitude, the Feynman diagrams of Fig. 121 must be 
calculated to O(p ), where p is the small momentum used as an expansion parameter. As 
discussed below Eq. (^), this implies that we need to include NNLO irnee operators, NLO 
NNnee operators and LO NNNNee operators. 

From the imee Lagrangian of Eq. (J3~2j) . the LO Oz//3/3-decay amplitude of Fig. |2k is 
calculated to be 

™ _ gjG 2 F A 2 H M 2 8(u p3 -f 5 u nl )(u p g 5 u n2 ) , 5 T 

M A 71 2 i ■ \f 2 o , - x X M el7 7 [Pi + P2l )U e2 , (.48 J 

A /3/? [Qi + ie) {qZ - mi + ie) 

where q\ = Pi — P 3 , q 2 = P 2 — P^ as defined in Fig. |2K and g& — 1.27 is the usual axial 
pion-nucleon coupling related to g^NN by the Goldberger-Treiman relation. 

As for the NLO, recall from Eq. (|15|) and the discussion that followed that the 7T7ree 
vertex has no NLO contributions. Thus, the NLO O^/3/3-decay nuclear operators are given 
by Fig. |2) andEt. Note that experiments planned and under way involve mainly ground state 
to ground state transitions + — > + which are favored by phase space considerations. The 
nuclear matrix elements of all the operators of C^ NlTee [Eq. (|39)) ] vanish for this transition 
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by parity 9 . There are therefore no NLO contributions for the + — > + transition and Mq 71 " 
is the only non- vanishing amplitude through 0(p). Nevertheless, we provide the expressions 
for the NLO nuclear operators in Appendix O for completeness. 

Taking the non-relativistic limit of Eq. (J48|) and Fourier transforming to co-ordinate space 
yields 

F.T.M ™ - T^^T^ n el 7 2 7 (A + ^7 5 )^ 2 ^(^^2,f3^4), (49) 

127T App 

where the nuclear operator is given by 

0^(X!, . . . , f 4 ) = -5(X! - X 3 )5(x 2 - ^4)(X3, Q Xl,/3)(xi,0X2,5)^[Fi(J a/3 ■ <t # + F 2 T ai>j p S ], (50) 

and 

T a <j>,p& = 3a Q/ g • pa^s • P — <J a p • &4>5 ■ (51) 

The form-factors Fi and F 2 were first introduced in Ref. j^J 

Fi(x) = (x - 2)e- x , F 2 (x) = (x + l)e~ x , (52) 

where x = m n p, p = \xi — x 2 \ is the distance between the nucleons, and p = pj p. However, 
in Ref. j^j], these form-factors were derived within a minimal extension of the standard 
model with only left-handed currents and heavy Majorana neutrinos; as was shown above 
by considering the possible representations to which the product of two left-handed weak 
currents can belong, this minimal extension cannot give rise to the LO niree vertex that 
yields these form-factors. In contrast, the derivation of F\ and F 2 was performed here by 
considering the symmetry properties of the quark operators that could generate the hadronic 
0i//3/9-decay operators without specifying the short-distance physics responsible for Qv(3(3- 
decay. 

Up to NLO, the O^/3/3-decay half-life is therefore 

— ~ — ^^p> I dE^Z + 2, E X )F{Z + 2, E 2 ) 

2 i o2\ m TTi „ 771 t al o2\ 2l| \ a \2 



21/2 1447r 5 ln2^ A m JWn 

1 



2 m + ti)PiE lP2 E 2 - (fit - P 2 2 ) PlP2 K]\M \ 2 , (53) 



9 Recall from above that N^ 5 N and Nj l N, i — 1,2, 3, are NNLO operators that couple the large and small 
components of the nucleon spinors. 
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where F(Z, E) is the usual Fermi function describing the Coulomb effect on the outgoing 
electrons with 

M = < V A ,z+2\ Y] — [Fi(xij)ai ■ Bj + Fix^T^r^ A ,z >, (54) 

Tij = 3<7i • pij&j ■ — 3i ■ (jj, (55) 
E 2 = E pp - E 1 , Pl = yjEf-ml . (56) 

Here pij is the distance between the i'th and j'th neutrons in the initial nucleus \^a,z > 
or the distance between two protons in the final state |^a,z+2 >, m e is the mass of the 
electron, R is a scale taken to be of the order of the nuclear radius 10 &iu\ acts on the spin of 
the i(j)'th neutron and the isospin matrix turns the i(j)'th neutron into a proton. Note 
that independently of the nuclear matrix element, the j3f — /?f part of the rate in Eq. (JH3J) 
is always considerably smaller (by at least a factor of ~ 10 from the kinematics) than the 
j3\ + /3| part which is the only one usually considered. 



A. Long Range Operators At NNLO 

Consider now the long range operators at NNLO. We are interested in comparing the LO 
and NNLO tree-level long range contributions and for simplicity we will ignore contributions 
from loops, m 2 insertions and the four-nucleon vertex which also contribute at NNLO 11 . 
Thus, we only need the hadronic operators of Eqs. l)H5j) and ({Hj) rewritten here 

M 2 = p- {f^-d^n-eiPs + /? 4 7 5 )e c + / ff A H ^ 5 r + 7r-JVe(C 5 + Ce7 5 )e c + h.c.} . (57) 

A/3/3 

The diagrams of Fig. |2k,b and c can be evaluated using the operators of Eq. (|57|) . The 
Fourier transform of the final result is: 

M2 = ir&ir- fci-rVCft + ^W e2 ortfu . . . , f 4 ) 

+« e i7 2 7°(C 5 + C67 5 )< 2 0^ NN (x u . . . , X A )] , (58) 

10 This scale is inserted to make the operator in Eq. (|54|l dimensionlcss. It is canceled by a corresponding 
factor of 1/R 2 in the rate. 

11 We also ignore recoil order corrections from the amplitude of Fig.|3i where is of O(p ). In this case, 
the rate will be dominated by terms in Eq. (|53[) . 
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with 



-S{x! - x 3 )S{x 2 - Xi){xl, a Xi,p){x\ y <pX2,5)-^ (59) 
~~J^ S ( Sl ~ ^(x 2 - xd(xt >a Xi,p)(xl,,pX2,s)-^ (60) 



and (a; = w^p as before) 



Q n7T 


-y(4-*)e-*, 






Lr 2 — 




2 + 2x + -x 2 - 
3 


-x 6 
3 


e~ x 


qttNN _ 


1 


-x 2 e~ x , 






qttNN _ 


-( 


l + x + \x 2 )e- x 







(61) 
(62) 

(63) 
(64) 

The new form-factors G^ and G™ stem from the iriree vertex while G^ N and G 2 NN (also 
given in Ref. [3|) stem from the NNnee vertex. In contrast to the zero-derivative case, the 
amplitudes stemming from these two vertices are of the same order in this minimal extension 
of the standard model. 

The corresponding half-life, assuming that Eq. (J58j) represents the only decay amplitude, 

is: 

1 1 fhc\ 6 gi G 4 Y 



T 



1/2 



64tt 5 ln2 V R J h A^c 4 



x 



dE^Z + 2,E 1 )F(Z + 2,E 2 ) 



7T7T 

2 



v^Ah 
9aM 



CM 



nNN 
2 



7T7V 

2 



v^Ah 
9aM 



( 6 M1 NN 



9aM 



9aM 



p 1 E 1 p 2 E 2 



Pip 2 m e 



(65) 



with 



M7 {wNN) 



<^ +2 iE(^) 3 [ G r ( ^ 7V) ( 



XijJfJi ' (Tj 



n -Kir(-KNN)( 



x^Ttj T+r+\q AZ >. (66) 
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We can compare the rates of Eq. (|53|) and Eq. (|65|) by assuming that all dimensionless 
constants are of the order of unity with 1 / ~ and Ah ~ 1 GeV, and that the nuclear 
matrix elements cancel in the ratio: 

Eq.© A 4 H , 

^ ZZ 1 « 10 3 . 67 

Note that this ratio agrees with our expectation based on power counting. We end this 
subsection by emphasizing that Eq. (|fi5|) is not the general formula for the 0z//3/9-decay 
half-life at NNLO (which must include all contributing terms including loops, recoil effects, 
NNNNee terms and corrections) since the LO contributions should be added if they 
do not vanish from symmetry considerations before squaring the amplitude. 

IV. PARTICLE PHYSICS MODELS 

While our discussion so far has been quite general and independent of the underlying 
physics of the lepton-number violation, we apply in this section our EFT analysis to two 
particle physics models: RPV SUSY and the left-right symmetric (LRS) model. 

A. RPV SUSY 

R-parity-violating supersymmetry can contribute to Ou (3 /3-decay through diagrams like 
the one in Fig. ^p. Since supersymmetric particles are heavy, their internal lines can be 
shrunk to a point in tree level diagrams yielding operators that involve only quarks and 
leptons. When the RPV superpotential is expanded to yield a lepton number violating 
Lagrangian, and a Fierz transformation is used to separate leptonic from quark currents, 



the result is 



21] 



Q2 

^T]q J T J T^ , (68) 



where 



Jp = qi°r + q, J s = qr + q, = qa^(l + 7 5 )r + g, (69) 



and r]g,r]j are quadratic functions of the RPV SUSY parameter, A' m defined in Ref. (21 1 



2tt \x in rM 

~9 G\m\ 



2a s — 

m- n 



21 



with k = qj. (70) 



Here M is the nucleon mass, rriq is a first generation squark mass, m~ g is the gluino mass, 
a s is the running SU(3)c coupling, and the + ■ ■ • indicate contributions involving the first 
generation sleptons and lightest neutralino 12 . Note that the dependence on Gp and M 
cancel from Eq. (|68|) . so that the effective lepton-quark O^/3/3-decay operator depends on 
five inverse powers of SUSY masses. 

It is useful to rewrite Eq. (|55|) in terms of our operators Of±: 



r e = ^ie(l + 7 5 )e c 



(71) 



The first thing to note is that C?Jj~ can be neglected for + — > + nuclear transitions. 
Secondly, from Table ITT1 we see that 0%+ gives rise to LO iriree and NLO NNiree operators 
and therefore contributes to the long range O^/3/3-decay operator of Fig. that is enhanced 
relative to the short range interaction of Fig. 1211 as observed by direct calculation in Ref . (21) , 
but derived with different assumptions about the scaling of the LEC. 

From Eqs. (J3Uj) and (|32j) . it follows that the LO iriree operator contributes dom- 

inantly to the Ovflfl-dec&y in RPV SUSY. The corresponding half-life formula is Eq. ([53)1 
with f3\ = $2 and with the substitution 

1 1 fA 



A m 4M V 7* + '"J- (?2) 

Obviously, a lower limit on the half-life can be interpreted as an upper limit on the coupling 
constants r]g and rjj. Making further assumptions about masses of SUSY particles, one can 
ultimately obtain model-dependent upper limits on the coupling constant A' m as discussed 
in Ref. 39| . 

Next, let us compare the scaling rules used here and in Refs. 21 1 and j^]. In the previous 



section, we used NDA to extract the relevant scales out of the dimensionful LEC's by using 



the scaling rule Eq. (|31J1 . The alternative method used in Ref. |2l| was to calculate the 
quark operator matrix element in the vacuum insertion approximation (VIA) and match 
the result to the hadron operator matrix element. 

Specifically, for the LO iriree operator of Eq. (|30|) we found that the dimensionful LEC 



12 The slepton/neutralino terms - which have complicated expressions - cause r/q ^ 77^. We have only shown 
the gluino contributions for illustrative purposes. 
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scaled as A H /^ while the VIA would predict 13 



LEC's ~ (n + \J P Jp\n-) « (tt + | J p |0>(0| J p |0 



4 

~ V \m u m +m d y ' (?3) 



where m u ^ are the light quark masses. Taking Ah = A x = 47rf n , the chiral symmetry 
breaking scale, and m u + = 11.6 MeV we find 

NDA {A^Uffl 



VIA 



0.7 . (74) 



J7T (m u +m d ) 2 

The NDA scaling is thus slightly smaller than that obtained from the VIA. Although they 



give results of the same order, VIA has proved to be unreliable in other contexts (see, e.g., 
the study of rare kaon decays in Ref. |40j) . We will therefore use NDA in what follows. 

Referring to Table |H1 it follows that there should be additional, subdominant contribu- 
tions from the operator iriree and from the NNiree operator at NNLO. The NNLO con- 
tributions from the NNiree vertex were considered in Ref. [3^ where detailed numerical 
evaluations showed that they contribute on average about thirty times less then the LO 
contribution. Our systematic analysis leads to the same qualitative conclusion (namely with 
regards to the NNLO suppression of p 2 / A H with respect to the LO), but differs from Ref. [3| 
in some respects. 

First of all, not all NNLO contributions were included. In particular, as pointed out 
above, the NNLO iriree operator contributes to O^/3/3-decay at the same order as the NNiree 
operator (called lir in Ref. 0) and the form-factors G\\ should be included. 

Secondly, our analysis shows that the NNNNee operator (the only one considered previ- 
ously in this type of analysis) gives contributions at NNLO. 14 In Refs. j^lj the suppression of 
that operator relative to the LO iriree contribution was only by a factor of ten for 76 Ge which 
is larger than what would be expected from our power counting (see also Ref. How- 
ever, this suppression is still in qualitative agreement with our analysis keeping in mind that 



13 Note that we do not take into account the color factor 8/3 of Ref. |21j since it is a number of 0(1) which 
does not involve any mass scale. It can therefore be absorbed in the LEC's which are undetermined. See 
also the footnote below Eq. Q. 

14 We note that the long-range operators considered in Ref. |4jJ through the induced pseudoscalar coupling 
terms of the nucleon current correspond to the NNLO contributions of Eq. 1)35(1 . The results presented 
by the authors of Ref. ^lj in the form-factor approach arc compatible with the EFT analysis given here 
since they only considered left-handed hadronic currents. 



23 



considerable uncertainty remains in the evaluation of nuclear matrix elements. Furthermore, 
although the traditional method of calculating the short-range Oz//3/3-decay operator using 
dipole form-factors may yield results of the correct order, the method is unsystematic 
with uncontrollable errors that cannot be easily estimated. 



B. Left-right symmetric model 



We consider LRS models that contain a heavy right-handed neutrino, and mixing between 
the right-handed and left-handed gauge bosons with g L m g R = g where g^ and #r are the 
left-handed and right-handed gauge couplings. The LRS Lagrangian is taken to be invariant 
under SU(2) L x SU(2) R x U(1) B _ L where B, L are the baryon, lepton numbers respectively. 
We will not be concerned with the CP-violating phases of the mixing matrix U R of the right- 
handed quark generations (the right-handed equivalent of the Cabbibo-Kobayashi-Maskawa 
matrix, denoted here U L ) nor the precise nature of the relationship between U R and U L (e.g., 
manifest versus pseudo-manifest LRS model) as the order of magnitude of the constraints 
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44 



45 



46]. 



obtained from experiments are broadly robust to the different possibilities 
We will use the standard Higgs sector composed of a left-handed triplet, Al, a right-handed 
triplet, Ar, and a multiplet, that respectively transform under SU(2) L x SU(2) R x U(1) B _ L 
according to (L, R,Y)= (3,1,2), (1,3,2) and (2,2,0). Their vacuum expectation values are 

/ n \ / o \ 





(Al) 








(Ar) 



V A L/ 



($) 



V A R y 




(75) 



Assume the following relation between the gauge and the mass eigenstates (ignoring the 
possibility of a CP-violating phase) 



W L = cosC^i + sinC^ 2 
W R = - sin (W ± + cos (W 2 , 
where ( is a small mixing angle between the mass eigenstates and, 

,2 



Mi 



iii 



9(2, /2 

2 1 



|(W.' 2 + 2A R 2 ), 



KK 

aT' 



(76) 

(77) 
(78) 
(79) 
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where Mw 12 are the masses of W\$- From these equations and the fact that 
\k 2 + k' 2 |/2 > \kk'\, we immediately obtain the important relation first derived in Ref. 47 1, 15 

Turning to experimental bounds on the masses and mixing angles, we will use for the lower 
limit on the right handed gauge boson M W2 > 715 GcV [48], which corresponds roughly to 



A < 10" 2 , (81) 



To put 
in Ref. 



imits on the mixing angle, we use recent results from superallowed + — » + /5-decay 
4fl| that imply a violation of the unitarity of the CKM matrix at the 95% confidence 
level. In the LRS model, unitarity can be restored by taking a positive value for the mixing 
angle with magnitude 

C = 0.0016 ±0.0007 , (82) 



V ud \ 2 + \V US \ 2 + \V ub \ 2 = 0.9968 ± 0.0014 , (83) 



given that one has 



in the Standard Model only |49| . A range of 2 x 10~ 4 <(<3x 10" 3 is allowed at 95% 
confidence level. Note that the discrepancy in the unitarity condition cannot be resolved 
by adjusting A because it enters the ordinary /5-decay amplitude quadratically and, thus, 
produces a correction smaller than 10 -4 [see Eq. (j81j) ]. In what follows, we will consider the 
range < ( < 3 x 10~ 3 and use the central value of Eq. ()82j) for some specific estimates. 
Note that for the central value of ( of Eq. (|82j) . we obtain an upper limit on M^ 2 from 
Eq. (jUJ) of 

M W2 < M Wl /y/C M w 2 < 2 TeV, for ( = 0.0016 . (84) 

With these bounds on Mw 2 and (, we can now estimate the relative order of magnitude of 
the graphs of Fig. HJ 

When the right-handed neutrino and Wl,r are integrated out, the amplitude of Fig. 0^ 
reduces to an operator of the form while Fig. |3Jd reduces to an operator of the form 
. 16 In previous treatments of O^/5/5-decay, only graphic with right-handed interacting 



15 From here on, £ will exclusively denote the magnitude of the mixing angle. 

16 Recall that the parity-odd LL/RR operator is suppressed at NNLO. 
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FIG. 4: Left-right symmetric model graphs. Fig. ^jp, involves the interaction of two right-handed (left- 



handed) currents while Fig.^fc depicts the interaction of left-handed and right-handed currents. 
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currents is considered and the impact of Wl-Wr mixing is neglected. Our analysis of the 
previous sections implies that the hadronic operators generated by 0%+ are suppressed by a 
factor of p 2 /Ay ~ 10~ 2 relative to those generated by Q\+ ■ Hence, taking into account the 
fact that the coupling of a (right) left-handed current with a (Wi)Wa involves a suppression 
factor of £ while a W2 internal line involves a suppression factor of A, we expect the tck 
operators generated by these quark operators to scale as 



(LR) 
) 



AC < 10" 



X 2 ^ < 10" 



(85) 



with all else assumed equal. Therefore, even if £ is ten times smaller than the central value 
in Eq. ()82jl . the contribution stemming from the mixing of left-handed and right-handed 
gauge bosons is still non-negligible. It may even be dominant. 



Such analysis may modify two constraints that relate the right-handed weak boson and 
neutrino masses, M W2 and M Nr respectively 17 . 

The first constraint stems from the requirement that the vacuum expectation value of 
the Higgs field Ar be a true minimum of the Higgs potential that generates the masses of 
the right-handed particles j^j]. The vacuum is then stable against collapse i 



stringent constraints on the one-loop corrections to the effective potential 
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nis imposes 
52 , bdl. In 



17 



For illustrative purposes, we assume the existence of only one right-handed neutrino. 
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FIG. 5: Constraints on the right-handed weak boson and neutrino masses (in TeV) in the LRS model. 
The solid lines stem from the vacuum stability (V.S.) constraint of Eq. H8ty) while the hyphenated lines 
correspond to limits imposed from Ov /3 (3-decay and Eq. i)8(J|) with the following values of mixing angle from 
longest to shortest dashes: Q = {3 x 10~ 3 , 1.6 x 10~ 3 , 0} with i = 1, 2, 3. Graphs (a), and (b) correspond to 
cases 1 and 2 of the text, respectively. Note that the value of the mixing angle (3 = cannot occur for case 
2 without simultaneously taking Mw 2 to infinity, while C2 corresponds to the central value obtained from 
CKM unitarity. The arrows indicate the lower bound M\\r 2 >715 GeV imposed by direct searches. The 
shaded, triangular regions in the graphs are the allowed values of the masses if the mixing angle is £2 • 
Mw 2 Mw 2 




M Nr 




particular, the loop corrections will involve terms of the form /cA^ln(A^/A^ 2 ) where k is a 
constant that depends on the particle masses. For the vacuum to be stable at large values of 
Ar, k must be positive to ensure that the minimum at the VEV is a true minimum and not 
simply a local minimum. The condition k > is equivalent to a condition on the masses. 
Following this formalism allows us to derive a relationship between My/ 2 and M^ R : 



1.65M W2 > M, 



(86) 



This constraint is represented in the graphs of Fig.^lby the fact that no value of (M Nn , M W2 ) 
below the solid lines is allowed 18 . 



A second relationship constraining Mw 2 and My p in t 
inferred from experimental limits on O^/3/3-decay [50, 



re LRS model with mixing can be 
from Eq. flSfl) with App = M Nn 



18 In Ref. I^H, the constraint that appears is 0.95Mw 2 > Mn^i the result of a typo p^ . 
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and choosing (3\ = ft 



2 



Q A/f 2 

l{>± ^ +{,)f< ^A^ ^' (87) 



l H' j 0k K vl I ^1/2 

G^ Z) = (G f cos^a) 4 (^ 



dExi^Z + 2, Ei)F[Z + 2, E 2 ) Pl E lP2 E 2 , (88) 



where defined by Eq. (|87j) . Ah ~ 1 GeV, Tjfy is the current limit on the half- 

life of the Oz//3/3-decay transition of a nucleus (A, Z) and where the functions Gq v were 
tabulated in Ref. 5(| for various nuclei. The matrix element M.\^ ,Z ^ is defined in Eq. (J54|) . 19 
In Eq. (JH7J) we have made explicit the scaling factors of Eq. (jHSJ) and also introduced a factor 

5 which parametrizes the p 2 /A H suppression of the NNLO 0^/3/3-decay operators relative to 
the LO operators. As mentioned above, the numerical evaluations in Ref. [3] suggest that 

6 ~ 1/30 which is the conservative number we will use. Thus, the A 2 term stems from the 
exchange of two W 2 s while the ( 2 term comes from the exchange of two W\s where (, being 
the magnitude of the mixing angle, is always positive. The relative sign between the (X and 
S(\ 2 + ( 2 ) terms on the LHS of Eq. (JH7J) cannot be predicted by EFT since we do not know 
the sign of the LEC's. 

For the values of half-life, G^' Z ^ and A4^' Z \ we will use the ones determined for 76 Ge 

Tyl > 1.9 x 10 25 yrs, (G^)' 1 = 4.09 x 10 25 eV 2 yrs, M$ c = 2, (89) 



where we extracted the value of Ai^ e from the value of A4 2n calculated in Ref. 1 39 1 and the 



limit on the half-life is at 90% confidence level . With these numbers, Eq. (jHZj) becomes: 

|CA ± 5(X 2 + C 2 )| < - Gc = ^ 10" 6 , (90) 



In the limit ( — > we obtain 

/ \ V 4 , X 1/4 

TeV. (91) 



/ . /3.8 TeV fi \ „, /TeV \ 



Our result is slightly smaller theu the resu.t obtaiued iu Refs. HQ for zero mixing aug.e. 
In Refs |5(| this constraint was calculated with the short range NNLO NNNNee operator 



19 From here on, we take cos 9q = 1. 
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of Fig. |21i using the dipole form factor approach. Note that we can reproduce exactly the 



values given in Refs. 
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by slightly adjusting the unknown constants /3i, in Eq. (|53j). 
To extract the constraint imposed by Eq. ( |9L)p on M Nr and M W2 , we need to consider 
three cases: 

• (1) the LO and NNLO terms have the same sign which corresponds to taking the plus 
sign in Eq. (|5Dj). 

• (2) they have opposite signs with (X > 5(\ 2 + ( 2 ), and 

• (3) they have opposite signs with (X < 5(X 2 + ( 2 ). 

We note that in all three cases, the upper limit on M W2 for ( > implied by Eq. (jSOj) always 
holds. 

Case 1: When solving the quadratic equation in A, we must keep the root that has the 
same limit as Eq. (J§U|) when 5, ( —> 0, 

C < A < ^ (-C + V(l - 45 2 )C 2 + 45i/) , (92) 

where we used Eq. (180)1 to obtain the first inequality. The first thing to note is that Eqs. 
(|90H92j) impose a lower-limit on the mass of the right-handed neutrino 

M^L > A /Mi 6(i + 26)( 2 = 1.8 , (93) 
TeV V 9 K ' 



assuming the central value of Eq. (J82J). This lower limit only depends on the mixing angle 
since 5 can in principle be calculated. In Fig. the constraint Eq. (}9"2*|) is plotted for three 
values of the mixing angle Q = {3 x 10~ 3 , 1.6 x 10~ 3 , 0}. 

In Fig. EK, we see that the larger the mixing angle, the larger the parameter space that 
is ruled out. In particular, for d, the largest angle that we are considering, the region 
allowed by Eqs. ()9*2*|) and ()9*3*j) is located below the constraint imposed by vacuum stability. 
Hence, a value of the mixing angle as large as d is excluded. In contrast, the central mixing 
angle value from CKM unitarity, (2, allows for a triangular region [bordered by the vacuum 
stability curve and Eqs. (|92j) ] of possible values for the masses. In particular, for £ 2 , we note 
that not only do we have the upper-limit of Eq. ()84)1 . but we also have Mw 2 > 1-6 TeV and 
Mn r < 3.2 TeV, which would constitute more stringent limits than that obtained from direct 
searches so far. For zero mixing angle, the entire region that is simultaneously above the 
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vacuum stability curve and the curve stemming from Eq. (|9*T|) is allowed. Thus, in general, as 
the mixing angle increases, the allowed region of parameter space shrinks while the minimum 
value of M\y 2 increases. The maximum mixing angle that results in a non-vanishing allowed 

2D 

region 

C < 2.2 x 10" 3 , with M W2 = 1.7 TeV, M Nr = 2.8 TeV. (94) 

Case 2: The condition of validity for this case, (\ > 5(\ 2 + ( 2 ), rules out the positive 
root of the quadratic equation in A, Eq. (J9T!|) . The limits on A are then 

C < A < i (C - V / Mf 7 4fo) • (95) 
We note that Eq. (|93j) imposes upper and lower limits on both and Mw 2 , 

10 6 (1 - 25) C 2 < |f < ^10^(1 - 45 2 )C 2 , < M W2 . (96) 

For £ 2 , we obtain in particular, 1.6 TeV< M^ R < 12 TeV and Mw 2 > 0.51 TeV. Note that 
the upper limit on M^ R for (2 is well above the constraint stemming from vacuum stability, 
Eq. (|8l)j). combined with the upper limit on M W2 given in Eq. (JH1J). Eqs. (jHSJ) also implies a 
new relationship between Mjy R and M^ 2 applicable only to case 2, 

where we neglected the 45 2 term. 

From the plot in Fig. the same analysis as in case 1 follows: as the mixing angle 
increases, the region of allowed values for the masses shrinks. As in case 1, (1 is already 
excluded while (2 allows for a triangular region of possible values for the masses. We note 
that Eq. ()97j) does not further constrain the allowed region of parameter space and has been 
included here for completeness. For this case, the maximum mixing angle is calculated to 
be, 

C<2.1xl0~ 3 , with M W2 = 1.8 TeV, = 2.9 TeV, (98) 

which are similar to the values found for case 1. 
20 Actually, a point in this case. 
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TABLE III: Order at which the left-right symmetric models with/without mixing and RPV SUSY 
contribute to the Of (3 /3-decay operators of Fig. |2 



Models 


Fig.Hi 


Fig.Eb,c 


Fig. Hi 


LRSM ( = 




P° 


p° 


LRSM C + 


p-2 


P° 


p o 


RPV SUSY 


p-2 


p-i 





Case 3: For the case (X < <5(A 2 + £ 2 ), we must keep the root that gives the correct 
upper-limit when ( — > since now the limit 5^0 cannot be taken. With the constraint 
on A stemming from the condition of validity of this case, (X < S(X 2 + ( 2 ), the inequalities 
satisfied by A are 

A < L (i - VT^W) , ^ (l + s/l^W) < A < 1 (C + v / (l-45 2 )C 2 + 4^) . (99) 

Thus, values of A located between the roots X± = C/(2<5)(1 ± yl — 45 2 ) are excluded. 21 
Note that for the two non-zero angles considered in Fig. the ranges defined by A > A + 
have already been ruled out by direct searches of right-handed bosons Q| and we are left 
with the first constraint of Eqs. (|9~H|) which does not depend on limits from 0z//3/?-decay. 
However case three appears to be entirely ruled out by Eq. (JHUj) . Indeed, approximating the 
remaining constraint of Eq. (|9*9*|) to A < we see that both constraints cannot be satisfied 
simultaneously. 

From Fig. Eland the three cases considered above, it follows that the effect of mixing on 
the mass constraint can be very important - a point not recognized previously. In particular, 
we see that non-zero mixing angles will generally exclude much of the parameter space by 
imposing much more stringent constraints on the masses and that the mass of the right- 
handed neutrino is bounded from below. We also note that quite generally, the mixing angle 
is constrained to be < 2.2 x 1CU 3 . 

We conclude this section by briefly comparing the left-right symmetric model and RPV 
SUSY. We observe that although both models can contribute to 0(p~ 2 ) to the operator of 
Fig. EK, only RPV SUSY contributes to Figs. Eb,c to 0{p~ l ) as discussed in the previous 

21 Since 1/(2(5) (( + s/Q- - 46 2 )C 2 + 4<5r/) > C/(25) (l - \f\ - 46 2 ) for all non-zero values of C and v, we 
need only be concerned with the A_ upper limit on A. 
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section. These results are summarized in Table ITTT1 

V. CONCLUSIONS 

Neutrinoless double beta-decay will continue to probe "new" physics scenarios that violate 
lepton number for some time to come. The existence of such scenarios is intimately related to 
the nature of the neutrino, namely, whether or not it is a Majorana particle. If a significant 
signal for Ou/3/3 decay were to be observed, one would know that the neutrino is a Majorana 
particle. However, one would not know whether the rate is dominated by the exchange of 
a light Majorana neutrino or by some other L-violating process that is also responsible for 
generation of the Majorana mass. Such L-violating processes could involve mass scales (A^) 
well above the weak scale. Thus, it is important to study the implications of Oz//3/3-decay 
for such scenarios - a task which we have undertaken in the present paper. 

In doing so, we have applied the ideas of EFT, which is appropriate in this case because 
there is a clear distinction of scales: A^ ^> Ah ^> p. We wrote down all non-equivalent 
quark-lepton operators of dimension nine that contribute to Ov (3 /3-decay, and showed how 
to match them to hadron-lepton operators by using their transformation properties under 
parity and chiral SU(2). We then organized the hadron-lepton operators (nnee, NNnee 
and NNNNee) in powers of p/An and discussed how the symmetries determine the type 
of hadronic operators that can be generated by each quark operators. In particular, we 
demonstrated that the hadronic operators generated by the interaction of two left-handed 
or two right-handed quark currents are always of NNLO. We also showed that EFT can 
classify particle physics models of 0z//3/9-decay in terms of the hadron-lepton operators they 
can generate and to what order these operators enter. In particular, we found that left-right 
symmetric models with mixing can potentially and considerably modify existing constraints 
on the masses of the right-handed particles. Indeed, a non-zero mixing angle gives far more 
stringent constraints on the allowed values of the masses of right-handed particles including 
a correlation between the mass of the right-handed neutrino and the mixing angle. We 
also found that a necessary condition for the existence of a region of allowed values of M^ 2 
and Mn r is £ < 2.2 x 10 -3 . For RPV SUSY models, we have also confirmed the previous 
conclusion that the dominant contribution stems from the inree operator which leads to more 
severe constraints on the corresponding RPV SUSY parameters than traditionally believed. 
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More generally, with this EFT analysis and using Table ITTl it can be immediately known what 
hadron-lepton operators can be generated by any quark-lepton operators appearing in any 
particle physics model that gives rise to Oz^/3/3-decay, and to what order these hadron-lepton 
operators will contribute. Finally, we note that deriving detailed information about a given 
scenario for L- violation will require combining information from a variety of measurements. 
As our analysis of the left-right symmetric model shows, using studies of Oz//3/3-decay in 
conjunction with precision electroweak measurements (e.g., light quark /3-decay) and collider 
experiments can more severely constrain the particle physics parameter space than can any 
individual probe alone. Undertaking similar analysis for other new physics scenarios and 
other probes of L- violation constitutes an interesting problem for future study. 

We thank P. Bedaque, M. Butler, R. Mohapatra, and M. Savage for useful discussions. 
P.V. thanks Prof. J. Hofejsi for his hospitality at the Center for Particle and Nuclear 
Physics, Charles University, Prague, Czech Republic. This work was supported in part 
under Department of Energy contracts DE-FG02-00ER41146, DE-FG03-02ER41215, DE- 
FG03-88ER40397, and DE-FG03-00ER41132 and NSF award PHY-0071856. 



APPENDIX A: EQUIVALENT AND VANISHING QUARK OPERATORS 

All operators proportional to e c 7 M e and e c <j^ v e vanish identically by virtue of the fact 
that the electron fields are Grassmann variables. For example: 

/ U \T 2 

= ze T 7 2 7°7 / "e 
= -e c 7/,e 

= 0. (Al) 

Note also that 7 5 cr^ 1/ = 2ie tJlUal3 a a /3 implies that e c 7 5 <T Atl/ e also vanish identically. In Ref . jl7 1 . 
these operators were incorrectly included in their super-formula 22 . 

22 However, they neglected them in their final analysis because they worked in the s-wave approximation. 
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Other color singlet operators that could potentially contribute to O^/3/3-decay are 

OU = (^V,a)(?tr + ?L , 6 ) = lof^, (A2) 

0++ = (<&T + a^q L , a )(<f R T + a^q L , b ) ± (^r+^ R , A (^ T +^ R , 6 ) = yO++ (A3) 

8 + + = (9£r+^ ?R , )(^ r +^ ?Lj6 ) = 0, (A4) 

Off" = (SrV^a + ftrV^)^ (A5) 

0&" = («tV^ g -&tV^)(#t + 7 ,^ (A6) 

1 =F a 

where the latin indices denote color and terms that involve the product of color octet currents 
are ignored (see below). Using Fierz transformations and the following formula, 

1 1 8 

Sabred = ^5 ad S cb + 2 X adKb , ( A7 ) 

1=1 

it is easy to prove Eqs. (|A2jl - (|A6)l . Note that the second term on the right-hand side of 
Eq. (|A7|) represents the product of two color octet currents. This term does not contribute 
since the asymptotic states are colorless and a completeness relation involving only hadronic 
states can be inserted between the currents. We therefore neglect this contribution. 

Even though two Fierz-related operators can arise due to different short-distance dy- 



namics, they are physically indistinguishable. Note that in Ref. 
operators were included as separate operators. 



17} . these indistinguishable 



APPENDIX B: NAIVE DIMENSIONAL ANALYSIS SCALING RULE 

To determine the scaling rules of the various fields appearing in the chiral Lagrangian, 
start with the relation between the axial current and the pion decay constant 4p| . 

(0\A a ^\Ap)) = Wf^, (Bl) 

which implies that tt is normally normalized by f n . Recalling that chiral perturbation theory 
is an expansion in powers of p/A-a, we scale pion derivatives by A R noting that pion loop 
corrections will involve factors of p 2 /(4irf n ) 2 ] this suggests that A R ~ 47r/ T . 

Since the action is dimensionless, we also have from the kinetic energy term of the pion 
field 

/ d 4 *^vr • = / d 4 x(A H /.) 2 ^ • ^ • (B2) 
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This shows that we can associate with d 4 x the scale (A.ufn) 2 - This is the origin of the last 
factor of Eq. From the parity-conserving pion nucleon coupling, we have 



/ d 4 x^iV 7 5 jhcN = [ d 4 x(A H /.) 2 /^iV7 5 /-f N. 

J U J A HJn A Kjn 



(B3) 



This shows that we can associate the scale Ah/ 2 with NN. 

Next, we note that since the axial current at the quark level is given by q^^q while 
a contribution to the axial current at the hadronic level is N^j^N, we can also associate 
with qq the scale Ah/ 2 . For a 0i//3/3-decay quark-lepton operator, this implies 

Gr 



1 dHqTq)(qT'q)(eT"e c ) = ^ I d 4 x(A H /^ 2 



(B4) 



AflpJ ™ ™ ' J v a '*' A H / 2 A h / 2 

Therefore, we can associate the scale G|/ 2 /A^g with the lepton bilinears. This explains the 
origin of the scaling rule in Eq. 1)311) . 



APPENDIX C: NLO NUCLEAR OPERATORS 



Here we present the results for Figs. |2b and I2t- The Lagrangian Eq. (J39j) gives 



(b) + (c) = At 



g A MA n _ 2 5 T 

—j=—u el -f 7 (Ci + C27 K2 



x 



(55 - m2 r + 

+ 4i^^M el 7 At7 2 7°7 5 Mj 2 x 



(gf - + z e ; 

(Mp37 5 M n l) 

(g 2 - + ie 



-u p3 (Cs + C47 5 ) + 
-M p4 (Ca + C4 7 5 ) 7 M M n2 



(CI) 



After taking the non-relativistic limit and performing a Fourier transform we obtain: 

F.T.dDU * f J^a 4^ *fo ~ ^)S(x 2 - <B,) — (1 + -) x 
2ttV2^ a Ah A^ p x 

M e i7 2 7°(Ci + C2 7 5 )mJ 2 (taS^ • pxi - Snxl° ■ PX2) 

+ n el7/ , 7 2 7 °7 5 'uj 2 



x 



(C2) 
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One can check explicitly that this nuclear operator is parity-odd and does not contribute to 
the + — > + nuclear transitions. Note also the extra factor of m n /An relative to the LO 
contribution of Eq. (|49|) which is consistent with the power counting of Eq. (JTJ. 



S. R. Elliott and P. Vogel, Ann. Rev. Nucl. Part. Sci. 52, 115 (2002) |arXiv:hep-p h/0202264 . 
J. D. Vergados, Phys. Rept. 361, 1 (2002). 

P. Vogel, Current aspects of Neutrino Physics, D.O. Caldwell (ed.), Springer, Berlin, 177 
(2001) |arXiv:nucl-th/0005020| . 

Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 81, 1562 (1998) 
|arXiv:hep-ex/9807003| . 

Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 87, 071301 (2001) 
|arXiv:nucl-ex/0106015| . 

K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 90, 021802 (2003) 
|arXiv:hep-ex/021 2021 . 

S. Pascoli and S. T. Petcov, Phys. Lett. B 544, 239 (2002) |arXiv:hep-ph /0205022| . 
R. N. Mohapatra and J. C. Pati, Phys. Rev. D 11, 2558 (1975). 

G. Senjanovic and R. N. Mohapatra, Phys. Rev. D 12, 1502 (1975). 
R. N. Mohapatra, Phys. Rev. D 34, 3457 (1986). 

J. D. Vergados, Phys. Lett. B 184, 55 (1987). 

M. Hirsch, H. V. Klapdor-Kleingrothaus and S. G. Kovalenko, Phys. Rev. D 53, 1329 (1996) 
|arXiv:hep-ph/9502385| . 

M. Hirsch, H. V. Klapdor-Kleingrothaus and S. G. Kovalenko, Phys. Rev. Lett. 75, 17 (1995). 
J. Schechter and J. W. Valle, Phys. Rev. D 25, 2951 (1982). 

H. Pas, M. Hirsch, S. G. Kovalenko and H. V. Klapdor-Kleingrothaus, Prog. Part. Nucl. Phys. 
40, 283 (1998) |arXiv:hep-ph/9712361| . 



H. Pas, M. Hirsch, H. V. Klapdor-Kleingrothaus and S. G. Kovalenko, Phys. Lett. B 453, 194 
(1999). 

H. Pas, M. Hirsch, H. V. Klapdor-Kleingrothaus and S. G. Kovalenko, Phys. Lett. B 498, 35 

(2001) |arXiv:hep-ph/0008182| . 

J. D. Vergados, Phys. Rev. C 24, 640 (1981). 



36 



[19] B. Pontecorvo, Phys. Lett. B 26, 630 (1968). 
[20] J. D. Vergados, Phys. Rev. D 25, 914 (1982). 

[21] A. Faessler, S. Kovalenko, F. Simkovic and J. Schwieger, Phys. Rev. Lett. 78, 183 (1997) 



|arXiv:hep-ph/9612357| . 
[22] U. van Kolck, Phys. Rev. C 49, 2932 (1994). 

[23] J. L. Friar, D. Huber and U. van Kolck, Phys. Rev. C 59, 53 (1999) |arXiv:nucl-tn79 809065 . 

[24] P. F. Bedaque and U. van Kolck, Ann. Rev. Nucl. Part. Sci. 22, 339 (2002). 

[25] S. R. Beane, P. F. Bedaque, M. J. Savage, and U. van Kolck, Nucl. Phys. A 700, 377 (2002). 

[26] J. Gasser and H. Leutwyler, Annals Phys. 158, 142 (1984). 

[27] J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985). 

[28] D. B. Kaplan and M. J. Savage, Nucl. Phys. A 556, 653 (1993) [Erratum-ibid. A 570, 833 

(1994 ERRAT,A580,679.1994)]. 
[29] A. Manohar and H. Georgi, Nucl. Phys. B 234, 189 (1984). 
[30] M. J. Savage, Phys. Rev. C 59, 2293 (1999) |arXiv:nucl-th/9811087| . 



[31] S. Weinberg, Phys. Lett. B 251, 288 (1990). 
[32] S. Weinberg, Nucl. Phys. B 363, 3 (1991). 

[33] D. B. Kaplan, M. J. Savage and M. B. Wise, Nucl. Phys. B 478, 629 (1996) 
|arXiv:nucl-th/960'5002l . 

[34] D. B. Kaplan, M. J. Savage and M. B. Wise, Phys. Lett. B 424, 390 (1998) 

|arXiv:nucl-th/9801034| . 
[35] T. D. Cohen and J. M. Hansen, |arXiv:nucl-th/9 908049 . 



[36] E. Jenkins and A. V. Manohar, Phys. Lett. B 255, 558 (1991). 
[37] T. Becher and H. Leutwyler, Eur. Phys. J. C 9, 643 (1999) jarXiv:hep-ph/9901384| . 
[38] D. Lehmann and G. Prezeau, Phys. Rev. D 65, 016001 (2002) |arXiv:hep-ph/0102161| . 
[39] A. Faessler, S. Kovalenko and F. Simkovic, Phys. Rev. D 58, 115004 (1998) 
|arXiv:hep-ph/9803253] . 

[40] J. F. Donoghue, E. Golowich and B. R. Holstein, Cambridge Monogr. Part. Phys. Nucl. Phys. 
Cosmol. 2, 1 (1992). 

[41] F. Simkovic, G. Pantis, J. D. Vergados and A. Faessler, Phys. Rev. C 60, 055502 (1999) 
|arXiv:hep-ph/9905509] . 

[42] A. Wodecki, W. A. Kaminski and F. Simkovic, Phys. Rev. D 60, 115007 (1999) 

37 



|arXiv:hep-ph/ 9902453 . 

[43] M. A. Beg, R. V. Budny, R. N. Mohapatra and A. Sirlin, Phys. Rev. Lett. 38, 1252 (1977) 

[Erratum-ibid. 39, 54 (1977)]. 
[44] S. Abachi et al. [DO Collaboration], Phys. Rev. Lett. 76, 3271 (1996) |arXiv:hep-ex/9512007| . 
[45] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 74, 2900 (1995). 
[46] P. Langacker and S. Uma Sankar, Phys. Rev. D 40, 1569 (1989). 
[47] E. Masso, Phys. Rev. Lett. 52, 1956 (1984). 

[48] K. Hagiwara et al. [Particle Data Group Collaboration], Phys. Rev. D 66, 010001 (2002). 

[49] J. C. Hardy and I. S. Towner, Eur. Phys. J. A 15, 223 (2002). 

[50] R. N. Mohapatra, Phys. Rev. D 34, 909 (1986). 

[51] S. R. Coleman and E. Weinberg, Phys. Rev. D 7, 1888 (1973). 

[52] S. Weinberg, Phys. Rev. Lett. 36, 294 (1976). 

[53] P. Q. Hung, Phys. Rev. Lett. 42, 873 (1979). 

[54] R. N. Mohapatra, private communication. 

[55] M. Hirsch, H. V. Klapdor-Kleingrothaus and O. Panella, Phys. Lett. B 374, 7 (1996) 

|arXiv:hep-ph/9602306] . 
[56] M. Doi, T. Kotani and E. Takasugi, Prog. Theor. Phys. Suppl. 83, 1 (1985). 
[57] H. V. Klapdor-Kleingrothaus et al, Eur. Phys. J. A 12, 147 (2001) |arXiv:hep-ph/0103062| . 



38 



